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Abstract 

We generalize the notion of quasi-local charges, introduced by P. Tod for Yang- 
Mills fields with unitary gauge groups, to non-Abelian gauge theories with arbitrary 
£jT)! gauge group, and calculate its small sphere and large sphere limits both at spatial 

and null infinity. We show that for semisimple gauge groups no reasonable definition 
, yield conserved total charges and Newman-Penrose (NP) type quantities at null 

infinity in generic, radiative configurations. The conditions of their conservation, 
■ both in terms of the field configurations and the structure of the gauge group, are 

. clarified. We also calculate the NP quantities for stationary, asymptotic solutions of 

the field equations with vanishing magnetic charges, and illustrate these by explicit 
CN ' solutions with various gauge groups. 

1 Introduction 

Just a couple of years after the introduction of the concept of (non-Abelian) gauge fields, 
in his classical paper Utiyama [1] already intended to interpret general relativity as a 
gauge theory, in which the Lorentz group emerges as the gauge group of transformations 
taking locally denned orthonormal bases to orthonormal ones. Then M0ller [2] formulated 
general relativity as a dynamical theory of frame fields, and in this new form Einstein's 
theory showed indeed remarkable resemblance to Yang-Mills theories. This similarity 
of the two theories made it possible to import techniques developed in one field to the 
other (see e.g. [3J), yielded a deeper understanding of both (see e.g. [4, 5J), or moti- 
vated the search for a proper gauge theoretic foundation of general relativity (see e.g. 
[6, 7J). In particular, the concept of charge density of non-Abelian gauge fields appeared 
to be similar to the (ill-defined) concept of the energy-momentum density of the grav- 
itational field: the vanishing of the gauge- covariant divergence of the (gauge-covariant) 
non-Abelian charge 4-current (e.g. of the fermionic matter) in itself does not yield any 
gauge-covariant notion of conserved charges. 
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To overcome this difficulty, Tod [8] introduced the notion of quasi-local charges in 
non-Abelian Yang-Mills theories, following the analogous twistorial definition of quasi- 
local energy-momentum and angular momentum of Penrose in general relativity. His 
construction is based on the use of the holomorphic and anti-holomorphic cross sections 
of a Hermitian vector bundle. Hence the original construction could give a well defined 
notion of charge in theories with (pseudo-)unitary gauge groups (or subgroups of them). 
However, the Lorentz group is not such a group, thus, to be able to see the closest analogy 
between general relativity and proper non-Abelian gauge theories as it could be possible, 
it would have to be considered the latter with the Lorentz (or its covering group, the 
iSX(2,C)) gauge group. Therefore, we need to generalize the notion of charge of Tod to 
arbitrary gauge groups, or at least to the groups above. 

In general relativity the energy-momentum and (relativistic) angular momentum take 
the form of 2-surface integrals, like the charge in Maxwell theory. However, in Einstein's 
theory there are ten additional quantities (and in the coupled Einstein-Maxwell system 16 
quantities), taking the form of 2-surface integrals on spherical cuts of future null infinity, 
which are absolutely conserved. These are the Newman-Penrose (or NP) quantities [9, 10J. 
Though their meaning is still not quite well understood, in the linear approximation they 
appear to characterize the radiation profile of the incoming radiation, while in stationary 
configurations they are connected with the multipole structure of the source. Thus the 
analogy between Einstein's theory and the non-Abelian gauge theories raises the question 
if the NP quantities can be introduced in Yang-Mills theories, and whether or not they are 
conserved. Though there is an argumentation that the analogous NP quantities cannot 
be conserved [HI H2], as far as we are aware no such mathematical analysis of the question 
has been given. 

Although the charges and the NP type quantities are conceptually different physical 
quantities, technically they are similar. Thus it seems natural to consider both. Therefore, 
the aim of the present paper is twofold: first, to generalize Tod's definition of quasi-local 
charges to arbitrary gauge group and investigate the resulting expression in the standard 
small sphere and large sphere limits. The second is to form the analog of the NP quantities 
and clarify the conditions of their conservation. 

Mostly to fix the notation, in section 2, we review the theoretical background and 
recall some basic facts that we use. Then, in section 3, we introduce the general notion 
of the self-dual/anti-self-dual charge integrals. We will see that the notion of electric 
and magnetic charges can be formulated and studied in terms of these more naturally 
and easier than the original charge integrals of the field strength. In subsection 3.2 we 
summarize the key points of Tod's construction for Hermitian vector bundles, and its 
generalization for arbitrary vector bundles and connections is suggested in subsection 
3.3, using either anti-holomorphic or holomorphic cross sections of the vector bundle. 
Here it is also shown that for globally trivializable bundles the modified construction also 
gives vanishing magnetic charges. 

In subsection 3.4. the limit of the quasi-local charge integrals are calculated in three 
limits. First, the charges are associated with small spheres near a point. We show that 
in the first two non-trivial orders the electric charge integral is independent of the actual 
construction, i.e. has some universal nature. The solution of the field equations near a 
point with the necessary accuracy is given in Appendix 15.11 Then we calculate the spatial 
infinity limit. We show that under the standard boundary conditions for the gauge fields 
this limit is finite and constant in time, i.e. these define conserved total charges. In the 
case of the vanishing total magnetic charges the total electric charges in the holomorphic 



2 



and in the anti-holomorphic constructions coincide. We also calculate the null infinity 
limit of the general charge integrals. We show that there is no reasonable definition for 
the total charges at null infinity that, for semisimple gauge groups, could be conserved 
in general radiative configurations. There can be charge conservation only in special 
(e.g. not radiative) field configurations or for gauge groups whose Lie algebra contains 
a non-trivial solvable ideal, provided the charges are introduced in the anti-holomorphic 
construction. The holomorphic construction appears naturally at past null infinity. 

Section 4. is devoted to the Newman-Penrose type quantities. First we show that 
for semisimple gauge groups no reasonable definition yields conserved NP type quantities 
in the generic, radiative configurations, and their conservation can be expected only 
in special configurations or for special gauge groups. We also calculate them in special 
stationary field configurations and find that they are built from the expansion coefficients 
of the field strength with dipole and higher moments. To be able to do these calculations 
we had to increase the accuracy of the asymptotic solutions of the Yang-Mills equations 
of Newman [3 J by two orders. These solutions are given in Appendix 15.21 

The signature of the spacetime metric is —2, and in connection with the spin weighted 
spherical harmonics we use the conventions of |13| . 

2 The theoretical background 

Let E(M) be a real or complex vector bundle over the spacetime manifold M with an n 
dimensional real or complex vector space ¥ A as the typical fiber, let {E^}, A = 1, n, 
be a locally defined frame field on some open U C M, i.e. at each point p e U the vectors 
Ei, E A span the fibre at p. (Thus the capital Latin indices A, B, ... are abstract 
indices referring to the internal fiber space, while the capital boldface Latin indices are 
concrete name indices taking numerical values. Small Latin indices are abstract spacetime 
indices.) The dual frame field will be denoted by {E A }, and hence E^E A = 5g and 
E^E B = #a hold. Clearly, such a frame field is not unique, and we allow the vector 
basis to be transformed as E^ i— > E^A A B , where A A b is a function on U taking its 
values in some subgroup G of GL(n,M) or GL(n,C). This yields the transformation 
E A i — y .E a Aa B of the dual frame field on U, where A A bAa C = <5g. Thus we think of the 
set of the frames above as a principal fibre bundle P(M, G) with structure group G, and 
the matrices A a b provide the free action of G on P. 

Let a connection D e be given on E(M). This can be specified by its action on the 
vectors E^ of the frame field, yielding the connection 1-form := E A (D e E^) on U. 
The connection is a G-connection if A^, as a locally defined matrix valued 1-form, takes 
its values in the Lie algebra g of the Lie group G. The curvature (or field strength), also 
in the local frame field on U, is the g valued 2-form F A Bcd = VcA^ — V^A^ + A A q A^ B — 
A dc A % which satisfies the Bianchi identity (V d F A Bef +A A :: F c Bef -F A cefA^ B )5'^ f c = 0. 
F A Bcd '■= E^E^F A Bcd is already gauge-covariant, and globally well defined. 

The Yang-Mills equations with gauge group G and gauge-covariant source current 
j A Be are 

g cb {V c F A Bba + A A c F c Bba - F A c6aA c c B ) = -47Tj A B a. (2.1) 

These equations imply, in particular, that the spacetime and gauge-covariant divergence 
of j A Ba is vanishing, and hence its spacetime divergence in itself is not: V a j A Ba — 
gob^jA^^c^ _ A A c j c Ba ). Therefore, by taking the flux integral on some spacelike 
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hypersurface E, it is the V e -divergence free and not gauge-covariant current J A Ba '■= 
j A Ba + 4~(-^ A Cab^4dB — A A c F c Bab )g hc that should be used to define charge. Thus it is 
more convenient to use the non-gauge-covariant form 

V a F A B a6 = -4?r J A Bb (2.2) 

of the field equations. In the rest of the paper we frequently use the matrix notation 
instead of writing out the boldface (matrix Lie algebra) or abstract fiber indices explicitly. 
For example, equation (12. ip takes the form g cb (V c F ba + [A c , F^]) = — 4irj a . 

If we fix a basis {e^ B } in the Lie algebra g, a = l,...,dimg, then we can write 
A^q = A®e A B and F A Bc d = F a c de A B . We define the structure constant c^o of g in 
this basis in the usual way by e A B e B c — e A B e B c =: c^^e A c . Then the definitions yield 
F^ab — V ' a A^ — V b A^ + cf^oA^Afr, and that the gauge-covariant derivative of any field 
e.g. with components X A B ='■ X a e A B can also be written in the form D e X A B = 
(VeX^ + c^ i/3 A"X l3 )e A B . The expression between the parentheses is just the gauge- 
covariant derivative on the vector bundle defined via the adjoint representation of G 
as an abstract Lie group, i.e. in which the typical fiber is g as its abstract Lie algebra. 
Then it is straightforward to rewrite the Yang-Mills equations in terms of these notions. 
Also for later use, we introduce the symmetric metric G a p '■= e^ B e^ A on the Lie algebra 
g. It is a direct calculation to check that e A B (which are constant functions on the domain 
U of the frame field {Ej[}) as a field on U is constant with respect to the gauge-covariant 
derivative D e , too. Therefore, the connection is compatible with this metric: D e G a p = 0. 
In the adjoint representation, in which the basis {e A B } is chosen to be the structure con- 
stants d^g, this metric reduces to the Killing-Cartan metric g a p on g. g a p is well known 
to be non-singular precisely for semisimple g, and this is negative definite precisely for 
Lie algebras of compact semisimple (and hence real) Lie groups. 



3 Quasi-local charges 

3.1 The concept of quasi-local charge 

Let S be any smooth, orientable, closed spacelike 2-surface. Let us fix a Newman-Penrose 
complex null tetrad {l a , n a , m a , rh a } on S (i.e. l a and n a are future pointing real null 
vectors orthogonal to S, m a and rh a are complex null tangents of S, and they satisfy the 
normalization conditions l a n a = 1 and m a m a = — 1), and recall that the area 2-form on 
S can also be given as e a b = — 2im[ a fhb], an d the area element is dS = \s a b ( see - e.g. 
|14|). Note that while l a and n a are globally defined on S, m a and rh a are defined only 
on the local trivialization domains of the tangent bundle of S. Then we search for the 
electric and magnetic charges in the form 



£ H : = 4^ j s sA b * F B Aab m a m b dS = i- j Tr(e( X i + Xi))d<S, (3.1) 
M M := ^ jf fi A B F B Aab m a m b dS = ± £ Tr(fi( Xl - Xi))d5, (3.2) 

where e A b and h A b are two appropriately chosen globally defined cross sections of E(S) <g) 
E*(S), where E*(S) is the dual bundle, *F ab := \e ab cd F cd is the dual of the field strength, 
where e abc d is the spacetime volume 4-form. In the expressions on the right we adopted 
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the matrix notation, in which Xi '■= \Fab{l a n b + fn a m b ) is just the m a m b component of 
the ant i- self- dual part of the curvature, and Xi is the m a m b component of the self-dual 
part of the curvature (see e.g. [2J US]). (N.B.: *F a bfh a m b = —iF a bl a n b . Moreover, 
though Xi is defined as the rh a m b component of F a ^ and the vectors m a and rh a are only 
locally defined, Xi lS globally well defined. In fact, it can be rewritten into the form 
Xi = — jF a b(e + ie ab ), where Eab '■= n a lb — l a nb is the area 2-form on the timelike 
2-planes orthogonal to S. In terms of the electric and magnetic field strengths defined 
with respect to some timelike unit vector field t a orthogonal to S, E a := F a i,t b and 
B a := *F a f J t b , respectively, Xi = \{E a + iB a )v a holds, where v a is the outward pointing 
unit spacelike normal of S such that t a v a = 0.) For real bundle and connection Xi and 
Xi are complex conjugate of each other, but in the complex case these are independent. 
Note that for e A b = fJ> b = S B the first integral can be rewritten into the flux integral of 
the trace (in the Lie algebra indices) of the current J a on any spacelike hypersurface E 
whose boundary is S, while the second is precisely the integral of the first Chern class^ 
of E(S), the pull back of E(M) to S. Thus S[e] and -M[//] are interpreted as the electric 
and magnetic charges, respectively, surrounded by S and defined with respect to the fields 
s A b and \x a b- 

If we define q A B '■= £ A b — ^ A b, then 

Q[q] :=£[q]-iM[q] = ^ jf T%%i)d<S, (3.3) 

in which only the anti- self- dual curvature component Xi appears. Xi appears in Q[q] := 
£[q] + LM[(7] with q A B '■= £ A b + ifJ- A B- If the bundle E(M) and the connection D e are 
real and e A b and /j, A b are chosen to be real, then Q[q] is simply the complex conjugate 
of Q[q\. For complex bundles and connections Q[q\ and Q[q] are independent. As we will 
see, it is more convenient to work with these anti-self-dual and self-dual charges than the 
original £[e] and .M[/i]. 

Since s A b and fi A g are still arbitrary, the number of the 'charges' £[s] and A^[/^] is 
infinite. However, we expect that the number of the independent charges be the number 
of the gauge potentials, i.e. dimg. Thus e A b and fi A B should be restricted so that the 
number of independent charges be the correct ones, and the resulting integrals define 
reasonable electric and magnetic charge multiplets. 

3.2 Tod's suggestion 

Tod's suggestion |S] for the quasi-local charges has the form (I3.ip - (l3.2p with e A b = \i A b — 
4> A <f>B f° r appropriately chosen cross sections (j) A of E(S) and yr B of E*(S): suppose that 
E(M) is a complex vector bundle endowed with a non-degenerate Hermitian fiber metric 
Hab 1 , and assume that the connection D e on E(M) is compatible with this fiber metric, 
D c Hab' — 0. Thus the gauge group is assumed to be some (pseudo) unitary group (or 
a subgroup of such). Then by means of Hab 1 the complex conjugate bundle E(M) can 
be identified with E*(M) via the fiber map tp A ' h-> tp* B := Hba iL P A ' ■ Finally, the section 
(p A was considered to be anti-holomorphic and tp A to be holomorphic, i.e. on the local 
trivialization domains to be satisfying 

1 Since we adopt the convention for the wedge product in which the wedge product of the two 1-forms 
a a and j3 a is (aA[3) a b = \{a a Pb~ Pa a b) (rather than only a a fib — /3 a otb), the proper normalization factor 
of the £:th Chern class is -^(j^) k (instead of y {^) k )- 
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= m e D e( p A = E A (50 A + m e A A B ) , 
= m e D e tp A = E A (5V A + m e A A B tp B ) . 



(3.4) 
(3.5) 



Here 3 and 8' are the standard edth operators on S and (fi A and Lp A are of zero-spin- 
weight functions on the domain of the frame field (see e.g. [HJ HI])- If V is the space 
of solutions 4> A of (l3.4j) and U is the space of solutions <p A of (13. 5p . then Tod defines 
the quasi-local charges to be the eigenvalues of the bilinear mappings £ : V x U — > C : 
(cf) A ,ip B ) t-> S[(f) A (p B ] and M : V x U -)■ C : (4> A ,Lp B ) h-> .M^^Vjg]. In the generic case 
on S ~ S* 2 this construction gives n (not necessarily real) charges, independently of the 
dimension of the gauge group. The magnetic charges are identically vanishing for globally 
trivializable bundle E(S). (For the details see [8].) 

3.3 The modified construction 

Our aim is to modify Tod's construction in such a way that (1) the gauge group can 
be arbitrary, (2) the number of the electric and the magnetic charges (at least in the 
generic case) be just dimg, and (3) the charges be real for real bundles and connections; 
but at the same time keeping the advantageous properties of the original construction 
(especially the vanishing of the magnetic charges for globally trivializable bundles). 

In principle there are several ways to ensure the correct number of charges. The first 
is to consider those cross sections whose components e A b and /i A B in a local frame field 
{E^} (as real or complex n x n matrix valued functions on the domain of the frame field) 
take their values in the matrix Lie algebra g. In this case, we could always write them 
as £ a e A B and /i a e A B for some, still arbitrary (real or complex) valued functions e a and 
fi a . Then it would be the (all together 2 dimg) functions e a and fi a that would have 
to be restricted appropriately to depend only on 2 dimg parameters in the generic case. 
Though conceptually this appears to be the most natural approach, it turns out that this 
framework is less flexible as it is more difficult to prove statements (e.g. the vanishing of 
the magnetic charges for globally trivializable bundles) than e.g. in the third one below. 

Another strategy is to build the charge integrals from the curvature and the fields 
e A b and /j A b in the adjoint representation, i.e. when the typical fiber ¥ A of the bundle 
E(M) is g as an abstract Lie algebra, and it is the structure constants that provide the 
matrix representation of a basis of g. However, the kernel of the action of the Lie algebra 
on itself via the adjoint map is not trivial, it is the center of the Lie algebra. Thus this 
construction may work e.g. for semisimple Lie algebras, but it would not yield any charge 
e.g. for Abelian Lie algebras. 

The third possibility is to construct the charge integrals in a vector bundle based 
on a representation in which the dimension n of the fiber space is just dimg (though 
the representation is not a priori the adjoint one). Since typically the dimension of the 
classical Lie algebras grows with the square of the dimension of their defining representa- 
tion, the representation space with the dimension of the Lie algebra appears to provide 
enough room for a faithful representation. Restricting the original construction of Tod in 
this way we get the correct number of charges. For the sake of simplicity we follow this 
strategy, but, instead of e A b and (j, A b, we formulate our conditions in terms of q A B and 

We search for q A B and q A s m the form 4> a ^>b and <f) A if)B, respectively, where <p A and ifis 
are anti-holomorphic while 4> A and ifis are holomorphic. If the bundle and the connection 
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are real, then <p A and tpB are anti-holomorphic cross sections of the complexified vector 
bundles E(S)(g>C and E*(S)<S>C, respectively, and the tilde denotes complex conjugation. 
In terms of their locally defined components these conditions are 

50 A + m e A A B <j) B = 0, 3^a - m e Af A ^j B = 0, (3.6) 
5> A + m e Af B 4> B = 0, 5> A - m e Af A ^ B = 0. (3.7) 

Then clearly q A B = 4> a ^b is anti-holomorphic and q A s = (p A ifjB is holomorphic. 

Next let us calculate the integrand of A^[/i]. Since 2i/r 4 # = q A B — q A B, by forming 
total derivatives and using (13. 6p and ( 13. 7p . the integrand of -M[/t] is 

2i f i A B F B Acd m c m d = d'(Tr((q-q) 1 )^ - $(Tr((q - q)j)"j + Tr((3g)7 + (8'g) 7 ) , 

where 7 := m e A A s and 7 := rh e A A 3 , and, for the sake of brevity, we used the matrix 
notation. For general anti-holomorphic q A B and holomorphic q A B we could not prove 
more. If, however, we use the special form of q A B and q A s given in terms of <p A , ipg, (ft A 
and i/jb and using (13. 6p and ( 13. 7p again, we obtain 

2i// B F B Acd m c m d = i5 a (e ah 4> A (5 b ^ A - i^Af A ) - e ab <p A {5 b ^ A - ^4f A )), (3.8) 

where 5 a is the intrinsic Levi-Civita derivative operator on S. Thus if E(S) is globally 
trivializable and hence the connection 1-form A^q can be chosen to be globally defined on 
S, then the right hand side of (I3.8P is also globally defined, yielding identically vanishing 
magnetic charges. Therefore, the modified construction also satisfies one of the minimal 
requirement of reasonableness. 

Let % denote the space of the anti-holomorphic cross sections <p A of E(S), and T-L* the 
space of the anti-holomorphic cross sections ip A of the dual bundle E*(S). Similarly, let 
% be the space of the holomorphic cross sections (f) A and %* the space of the holomorphic 
cross sections if) a- Then the anti-self-dual and the self-dual charges are defined to be the 
eigenvalues of the C-bilinear maps 

Q : U x W C : (0 A , Vb) H> Q [<I> A M , (3.9) 

Q-.nxn* ^C: (0 A ,^ B ) h. q[Wb], (3.10) 

respectively. If E(M) and _D e are real, then the complex conjugate of an anti-holomorphic 
cross section is holomorphic, and hence H, is the complex conjugate of H. Similarly, %* 
is the complex conjugate of W. Therefore, Q is the complex conjugate of Q, yielding 
real electric and magnetic charges. Repeating the argumentation of Tod in |8J one can 
show that for globally trivializable bundles over S ~ S 2 in the generic case the number 
of the electric charges is dim g. On the other hand, there are exceptional 2-surfaces for 
which the construction still works, but the number of charges is different. 

Though in general q A & = 5 A is not an element of g, it can be considered as the 
generator of the central extension of the matrix Lie algebra g. Since it is both holomorphic 
and anti-holomorphic, it is worth considering this extension. £[S\ and Ai[S] can be 
thought of as the 'mean' electric and magnetic charge, respectively, and, as we already 
noted, the latter is just the first Chern class of E(S). 
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3.4 Three limits 



Though the quasi-local charge integrals Q[q] are well defined even in curved spacetime, 
for the sake of simplicity in the rest of the paper we assume that the spacetime is the 
Minkowski spacetime. In particular, we calculate their small and large sphere limits in 
the flat spacetime. 

3.4.1 Small spheres 

Let p G M be an arbitrary point, t a a future pointing timelike unit vector at p, and 
let us consider the future pointing null geodesies starting from p with afline parameter 
r and tangents l a satisfying l a t a = 1 at p. Then for sufficiently small r the set <S r of 
the points on these null geodesies whose affine distance from p is r is a smooth spacelike 
2-surface, which is homeomorphic to S 2 . Such a surface is called a small sphere about 
p. Our aim is to calculate the charge integral Q[q] for small spheres of radius r, and to 
determine the coefficients in its expansion Q[q] = + ■ • • + r 4 Q^ + 0(r 5 ). (Since 
S r is contractible to p, any bundle over S r is globally trivializable, and hence -M[/t] is 
identically vanishing.) The strategy is to expand the curvature component Xi an d the 
solutions <p and ip of ( 13. 6 p (or of ( 13. 7p ) by powers of r. Since dS r = r 2 dSi, clearly = 
and Q^ 1 ' = 0, and to calculate Q[q] with r 4 accuracy we need to know xii 4> an d "0 with 
r 2 accuracy. 

Similarly, expanding the component 701/ := A a m a of the connection 1-form as a series 
of r, substituting this to (13.61) and taking into account that 5 = - Q, where 08 is the 
edth operator on the unit sphere (and is given explicitly in the appendix), we see that 
we need to know the expansion coefficients 7^, for I < 2. However, as the analysis of 
Appendix 15.11 shows, to determine the solution of the Yang-Mills equations on the light 
cone (more precisely, on the null boundary of the chronological future) of p for Xi with 
r 2 accuracy, formally 701', 710' and 711/ must be expanded with accuracy r 3 . The solution 
of the hypersurface equations fj5.ip ~ fj5.3p is given in Appendix 15. 1[ by means of which the 
equations (13. 6p themselves are 

ty<°> = 0, O 30 (1) = O, o 30 (2) = -\xT<P {0) ; (3.H) 

o3^ = 0, 3^ (1) = 0, 3^ 2) = \^xf- (3.12) 

Since and ip are multiplets of type (0,0) scalars, (f>(°\ <jP~\ tp^ and tp^ are constant 
on S. 

To solve the third of (13. lip and (I3.12p first let us recall that in terms of the spinor 
form of the field strength, F a b = <&ab£a'B' + &A'B'£ab, and the GHP spinor dyad {o A , t A } 
the curvature component x^ i s given by §abo a o b . (In the present subsection we use 
two-component spinors, and hence, in particular, the capital Latin indices A, B, ... in 
this subsection are abstract spacetime spinor indices. These should not be confused with 
the indices referring to the internal fiber space or the concrete name indices A, B, ... 
etc.) However, in terms of the Cartesian spinor dyad {0 A , I A } at p the GHP dyad is 

o A = -^-L (co A + i A ) , l a = ir _ j . (0 A - (I A ) , 

and hence it is straightforward to integrate (13. lip and ( I3.12p in the coordinates (C,C) : 
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4> {2) = ~YT7c ( C$0 + 2$1 " ^ 2 ) 0(O) + ^' (3 - 13) 
¥ 2) = y^^ 0) (c $ o + 2$! - c$ 2) + 4 2) - (3-14) 

Here $ : = ®abO a O b , $1 := § AB A I B , $ 2 := § A bI A I B ] and 0[, 2) and t/^ 2) are arbitrary 
constant sections, being the general solution of the homogeneous equation. Then, using 
the pattern of (I3.13p -( 13.14p . it is easy to write down the holomorphic solutions, too. 

Though for small spheres the solution spaces H, and H* are dimg dimensional, the 
space of the 0(r k ) accurate solutions is (k + 1) dimg dimensional. The presence of the 
'extra', spurious solutions is a consequence of the lack of a natural isomorphism between 
the solution spaces on different two-surfaces, and hence, in particular, with different radii. 
These yield some ambiguity in higher accurate approximations. (For a more detailed 
discussion of this issue in general relativity, see e.g. j 16 j . Also, in the calculation of 
certain integrals, we will use equation (2.5) of |16|.) 

Since ft ^ and ift(°> are constant, by (15. lip the 0(r 2 ) order term is vanishing. For 
the 0(r 3 ) order term we obtain 

Q (3) = ^/(^ (0) xS°V (1) +^ (0) xiV ) +^ X fU {0) )dS 1 = y^(^0 (O) ), (3.15) 

where the integration is taken on the unit sphere. Thus is the volume of the 3-ball 
of radius r times the timelike component of the gauge covariant current at p. This result 
coincides with that of Tod, too. In fact, let us observe that in only the zeroth 
order terms of the cross sections (p and ip matter (the others integrate to zero), which are 
constant both for the holomorphic and anti-holomorphic sections. 
In the fourth order we obtain 

Q (4) = ^ M '^^^4^(^ 0, e^ (0) ) + (3.16) 
+ f {t AA 't BB ' + t AB 't BA ') (VoiO (^ (0) e^ (0) ) + y t a f a (V (1) e^ ) + V (0 W«) • 

Though formally is built from the 0(r 2 ) accurate cross sections, their 0(r 2 ) order 
correction terms are integrate to zero. Hence we obtain the same result in the holomorphic 
and the anti-holomorphic constructions. On the other hand, in general does depend 
on the 0{r) order 'gauge solutions' ft 1 ' an d ifj^: is well defined only if = 0, 

i.e. when j£t a = at p. 

The calculation of is much longer and technically more difficult, so we summarize 
only the main message of that. First, ambiguities similar to that in parameterized 
by the 'extra' solutions ft 1 ', ... etc, appear, unless the third and fourth order terms, 
and are vanishing. The difference of the holomorphic and anti-holomorphic 

constructions can appear in this order. Thus and have some 'universal' nature, 
not being sensitive to the details of the defining equations for the cross sections and ip, 
provided <fi and ip are constant in the zeroth and first orders. 
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3.4.2 Large spheres at spatial infinity 

Let (t,r,(,C) be the standard spherical, complex stereographic coordinate system (at 
least in a neighbourhood of spatial infinity) in Minkowski spacetime, and consider the 
charge integral Q[q] on the metric 2-spheres <S r of radius r in some t = const spacelike 
hyperplane. Let t a and v a be their future pointing unit timelike and outward pointing unit 
spacelike normals, respectively, such that t a is orthogonal to the t = const hyperplanes 
and t a v a = 0. Let m a and rh a be the complex null tangents of the 2-spheres. Then we 
impose the following a priori fall-off conditions for the scalar potential and the radial and 
tangential parts of the spatial vector potential: $ := A a t a = 0(r~ a ), A a v a = 0(r~ b ), 
7 := A a m a = 0(r~ c ), 7 := A a rh a = 0(r~ c ), respectively, for some a,b,c > 0. In addition, 
we assume that the time derivative of the last three satisfy (A a v a ) = 0(r - ^), 7 = 0(r~ 7 ) 
and 7 = 0(r~ J ) also for some /3, 7 > 0. By technical reasons we also assume that the 
differentiation with respect to r reduces these orders by one, but the differentiation with 
respect to £ and £ does not change these orders. Then, expressing Xi — \ (E a + iB a )v a 
by the potentials, we find that Xi = 0(r~ 2 ) if a = b = c = 1 and (3 = 2. Therefore, these 
boundary conditions ensure the existence of the r — > 00 limit of the charge integrals 
Q[q] on 2-spheres S r of radius r provided q = 0(1). (If q were not involved in the 
2-surface integral of the curvature, i.e. if that were only the 2-surface integral of the 
components of the field strength in some globally defined frame {E^}, then one could 
rewrite that as the 3-volume integral of its divergence on some spacelike hypersurface E 
with boundary S. Then, using the constraint equation, one could show that even weaker 
boundary conditions could ensure the finiteness of the original 2-surface integral \17\ . 
Since, however, q is defined only on S, here this strategy cannot be followed, and hence 
we have the above 'conservative' fall-off conditions.) 

Next let us determine the conditions that ensure the vanishing of the r — y 00 limit of 
the time derivative of Q[q], 

j t Q[q]=—jTr[qXi + q^{E a + iB a )v a yS. (3.17) 

To determine the order of the second term of the integrand we use the evolution equations 
for the field strengths in the present geometrical situation (i.e. flat spacetime foliated by 
spacelike 3-planes with vanishing shift vector). These are 

E a = (p c B d )£ cd ab t b — [$, E a ] , B a = -(3 c E d )e cd ab t b -[<t>,B a ], 

where D a denotes the space- and spatial gauge-covariant derivative operator (i.e. for 
example its action on any cross section <p A = E^<p A of E(M) near the spatial infinity is 
B a (f> A = E A (d b (p A + A^0 B )P a b and P b := 8 b a - t b t a is the orthogonal projection to the 
t = const hyperplanes). By these equations Xi — 0(r~ 2 ~ e ) if 7 = 1 + e. Thus if e > 
and q = 0(r~ € ), then the r — > 00 limit of the time derivative of Q[q] is zero. Note that 
by 7 = 1 + e the asymptotic form of 7 (and of 7, too) is 

7=y ) (C,C) + ^ 7 7 (1) (U,C)+o(r- 1 -), (3.18) 

i.e. its leading order term does not depend on t. 

Finally, let us write q = + r~ e q^ + o(r~ e ) and substitute into the equation of the 
anti-holomorphicity of the cross section q. We find 
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= 5g+ [7,9] = 

= ;(W 0) + [7 (0) ,9 (0) ]) + ^(W 1} + [7 (0) ,g {1) ] + [7 (1) ,? (0) ]) +o(r^). 

Thus its solution is such that the leading term does not depend on time, and hence 
q = 0(r~ e ). Therefore, under the fall-off conditions for the vector potential above, the 
spatial infinity limit of Q[q] exists and is conserved in time. It might be interesting to 
note that here we used only the anti-holomorphicity of q A B, but did not assume that it 
has the dyadic product structure § a i\)b- 

Note that in general 7^°- ) cannot be transformed by gauge transformations to zero 
globally unless the leading term of the v a component of the magnetic field strength is 
vanishing. Thus q^ is not constant on the unit sphere, and hence the total charge depends 
on whether q is anti-holomorphic or holomorphic. If the connection coefficients 7 ( - ^ and 
7^ can be made vanishing globally, then the total magnetic charges are vanishing and the 
total electric charges introduced in the holomorphic and anti-holomorphic constructions 
coincide. 

3.4.3 Large spheres at null infinity 

Let us consider the charge integral Q[q] on any u = const cut of for still arbitrary 
globally defined q A s- (For a summary of the relevant results on the asymptotic structure 
of the Yang-Mills fields near J^ + , see Appendix 15.21 ) Because of the peeling this integral 
can be written in terms of the expansion coefficient Xi as j~ § Tr(?Xi)d<Si, where the 
integration is taken on the unit sphere. Hence this integral is still a function of the 
retarded time coordinate. In this subsection we show that for semisimple gauge groups 
no reasonable definition yields conserved total charge in general, radiative configurations. 
In particular, neither Tod's original definition (as pointed out in [S]) nor the modified 
one yield conserved total charges. 

By integration by parts and using (15 . 15 j) . for the ?/-derivative of Q[q] we obtain 

Q[q] =^j^(qx\~{^q+ [l° , q])xl)dS 1 . (3.19) 

Since in any given retarded time instant u = Uq in a generic (radiative) configuration Xi 
and X2 are independent data (because xi + X? an d 7° can be specified independently, see 
Appendix 15.2. ip , the integrand can be zero precisely when the coefficient of Xi an d of X2 
are vanishing. In particular, q A B must be anti-holomorphic, but in addition it would have 
to be constant in time even though the anti-holomorphic structure might be changing 
from cut to cut on 

q = 0, 3g+ [ 7 °,g] =0. (3.20) 

The condition of the integrability of the system of these equations is [q, 7 ] = 0. Thus the 
charge Q[q] with anti-holomorphic q A s is conserved in the absence of outgoing radiation, 
7 = 0, or when q is anti-holomorphic and belongs to the center of the Lie algebra g, too. 
Therefore, for semisimple gauge groups in the generic radiative configurations there is no 
total conserved charge. Charge conservation requires the non-triviality of the radical (i.e. 
the maximal solvable subalgebra) of g in its Levi-Malcev decomposition. 
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Remarkably enough, though in the presence of outgoing radiation there is no charge 
conservation, the structure of the Yang-Mills fields at future null infinity singles out just 
the anti-holomorphic cross sections q A s in a natural way Thus our suggestion for the 
quasi-local charges is justified by the conservation properties of the total charges at future 
null infinity At past null infinity the holomorphic cross sections would appear. 

Even for semisimple Lie algebras and non-stationary, but special configurations, in 
which Xo an d Xi + Xi are n °t quite independent, the total charge can be conserved. Such 
are the Lienard-Wiechert type solutions. In [18j Tod defined them geometrically as the 
algebraically special solutions of the Yang-Mills equations for which the field strength 
components Xo an d Xo are vanishing, and, in contrast to our present framework, the 
origin of the coordinate system (u,r,(,C) (i-e. the 'source's world line') is allowed to be 
an arbitrary timelike curve. He showed that there is an anti-holomorphic frame in which 
X? is constant both in time and on the cut of J^ + , and hence that the total charge is, 
in fact, conserved. (For another approach to the Lienard-Wiechert type configurations, 
based on certain assumptions on the explicit form of the 4-potential in terms of a pointlike 
source, the so-called colour, see [T^ [2D], and for a recent review of them see e.g. |21|). 

In the special stationary case discussed in Appendix 15.2.21 7 and 7 are vanishing 
and Xi — \C^e^B f° r some rea l or complex constants C M . Thus, in particular, the anti- 
holomorphicity of q A B yields that it is constant, i.e. q A B = Q^^aB f° r some constants q^. 
Therefore, Q[q] = C^e A B q B a = C a G a pq^ ', yielding the meaning of the constants C M of 
the solution: they should represent the total charge of the system. In fact, if we search for 
q A B in the form (j) A tpB for anti-holomorphic cf) A and ipB, then their anti-holomorphicity 
yields that they are constant, and hence that Q[q] = C u l e A B (j) E ^ a- Thus, indeed, the 
charges introduced in subsection 13.31 are the eigenvalues of C^e A B . 

4 Newman— Penrose type quantities 

The Newman-Penrose conserved quantities in general relativity and in Maxwell theory 
(and, in fact, in any theory based on a linear zero- rest-mass field equation with any 
spin or on the conformally invariant wave equation) are built from the second nontrivial 
expansion coefficient of the tetrad component of the curvature with the highest spin- 
weight, weighted with appropriate spin- weighted spherical harmonics [9j[T0]. Thus it is 
natural to define 



where a; is a cross section of (in the case of real bundles, the complexified) E(M)®E* (M) 
over the u = u cut of with —1 spin weight so that the integrand of 04. ip is globally 
well defined. 

4.1 The (non-) conservation of the Yang Mills NP quantities 

Here we show that for semisimple gauge group and any non-trivial choice for u the quan- 
tity Q[u] is not conserved in generic, radiative configurations. Using (I5.17P a straightfor- 
ward calculation (by integration by parts) yields 




(4.1) 
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Since Xo — Xo(CjC) an d Xo — Xo(C;C) are freely and independently specifiable functions 
at u = uq (see Appendix 15. 2.ip . Q[u] is constant in time precisely when their coefficients 
in (14. 2p are vanishing: 



uj = 0, 



a( 5W[7V])+ 7 ,(o5^+[7 ,c]) =3[ W>X ?]. 



(4.3) 



Using f 1 5 . 1 4 p - f 1 5 . 1 5 p of the appendix, the integrability condition of this system of equations 
is 



Thus, in particular, f2[u;] is conserved if a; satisfies (14. 3 j) and there is no outgoing radiation, 
i.e. %2 = 0; or ^ w belongs to the center of the Lie algebra g and solves o^o^'o; — 3[o>, Xi] — 
0. Therefore, the NP type quantities could be conserved in special situations (e.g. for 
a collection of independent Maxwell fields), but not conserved in general. However, this 
non-conservation of the NP quantities in non-Abelian Yang-Mills theories appears to 
be known for a while, as Bazanski pointed out to Newman in an unpublished private 
communication IT2]. 

Since for the Lienard-Wiechert type solutions of the Yang-Mills equations, considered 
in [T8], the curvature components Xo an d Xo are identically vanishing, the NP quantities 
for these solutions are vanishing for any gauge group. 



4.2 Yang— Mills NP quantities in stationary configurations 

It is known that in stationary configurations the NP quantities in the theories based 
on the linear zero-rest-mass field equations (e.g. in particular in Maxwell theory or the 
linearized Einstein theory) are vanishing. On the other hand, in stationary spacetime the 
gravitational NP quantities of the full non-linear theory reduce to a nontrivial combination 
of the mass and the dipole and quadrupole moments of the source [5J- Thus their non- 
vanishing is due to the nonlinearity of the theory. This result raises the question whether 
the analogous quantities Q[ui] in non-Abelian Yang-Mills theories could be no n- vanishing 
in the stationary configurations, and if they could, then what would be their meaning? 
In the special stationary case discussed in subsection 15.2.21 equation ( I5.18P reduces to 

2x1 = -oo^ - 5 [xi, Xo] " [o3Xo, Xo] - 2 [xo, o^Xo] ~ 3 [xi 0^X0] , 
while the defining equation (14.31) of uj to 

o a'a; = 3[a;,x?]. (4.4) 
Writing co = 0J^e^ B for some spin weight —1 functions and expanding the latter as 
= Yl™=x Em=-j u jm -i Y jm, equation flOD yields 

C\^ im := (C-C)^ = i (j - 1) (j + 2)u} m . (4.5) 

Thus can always be non-zero, and if C^ U) the 'charge matrix' in the adjoint repre- 
sentation, has only one eigenvector with zero eigenvalue then u^ m = C^uj m for arbitrary 
constants uj m . If, in addition to C M , there are other eigenvectors with zero eigenvalue, say 
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C^, ... etc, then u;f m will be a linear combination of all them with arbitrary coefficients 
(jj m , u m , ■■■ etc. However, u^ m for j > 2 can be nonzero only in the exceptional case 
when |(j — l)(j + 2) is an eigenvalue of CV (For a more detailed discussion of C^y see 
subsection 15.2.21 ) Substituting these into the definition of Q[u], by integration by parts 
we obtain 

Q[u] = ~ j u»G^c u a p(AC a kP + ( $h a )h p + 2h a ( Ni p ) + 3/i a ( 8V))dS 1 , (4.6) 

where the functions h^, h} 1 and k^ have been introduced via Xo =: ^ e «Bi Xo =: ^ e ^B 
and Xq =: k^e^ B , respectively. Thus, in particular, if the Lie algebra g is Abelian, then by 
( 14. 6 p Q[u] — 0, in accordance with the known fact that the NP quantities are vanishing 
for (an Abelian multiplet of) stationary Maxwell fields. Thus in the rest of the paper we 
assume that g is not Abelian. 

First suppose that C^ u is generic (see Appendix I5.2.3p . Then the only non-zero 
expansion coefficients of in terms of spin weighted spherical harmonics is u^ m , while 
the expansion coefficients of /i M and (and in the complex case of h^, too) are determined 
in Appendix 15.2.31 Substituting those solutions into equation ( 14. 6ft . using the notation 
of the Appendix and elementary properties of the edth operators and the spin weighted 
spherical harmonics (e.g. the orthogonality of these harmonics with different indices), we 
obtain 

1 r 1 

J m=-l 

In the special case when the only eigenvector of C^ v with zero eigenvalue is C M and — ~ 
is not an eigenvalue of C^ u this reduces to zero. Thus the non-Abelian nature of the 
Lie algebra g in itself is not enough to have non-trivial NP quantities. If, however, — ~ 
is an eigenvalue and k^ m are the corresponding eigenvectors, then Q[oj] can be non-zero. 
In the other extreme case when C^ v = (e.g. when the total charge is vanishing) it 
gives \G tlv ^2m= ) m ^fm w i-m) which is a homogeneous quadratic expression of the 
coefficients and characterizing the dipole structure of the Yang-Mills field. 

We saw in Appendix 15.2.31 that for the Lie algebras so(3), s/(2,M), and s!(2,C) the 
charge matrix has only one eigenvector with zero eigenvalue while for so(l,3) there are 
two, moreover there could be exceptional configurations for the last three Lie algebras, 
but not for so(3). Therefore, for the Lie algebra so(3) in general and for sl(2, R), sl(2, C) 
and so(l,3) in the generic case the NP quantities can be non-vanishing e.g. when the 
total charge is zero. 

Finally we consider the exceptional configurations of Appendix 15.2.31 In this case 
the solution of (@~5]) is = J2ln=-l U\ m -i*im + Y^L=-j u Jm -i Y Jm- Substituting this 
and the solution of (I5.33P into ( 14. 6p . and using the orthonormality of the spin weighted 
spherical harmonics on the unit sphere, we obtain that 

1 1 J 

m=—l m=~J 

which can also be non-zero in general, even if the charge matrix has a single eigenvector 
with zero eigenvalue. In this case, however, is built from higher, actually the 2 J th 
order multipole moments hj m and hj m of the curvature. 
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5 Appendix: The solution of the Yang-Mills equations 



In this appendix we summarize the structure of the Yang-Mills fields near a point (sub- 
section [5]T]) and near the future null infinity in Minkowski spacetime (subsection 15. 2p . For 
the generalization to the coupled Einstein- Yang-Mills equations in the Newman-Penrose 
formalism, see |15| . 

Following j3] we use the coordinates (u, r, CjC)' where u := t — r is the standard 
retarded null coordinate and (C, Q are the complex stereographic coordinates on the 
unit sphere. The complex null Newman-Penrose tetrad {l a , n a , m a , rh a } is adapted to 
the u = const, r = const 2-spheres, the edth operators on the unit metric sphere will 
be denoted by and qQ' and their action on a scalar rj of GHP type (p, q) is := 
^(1 + COg + i^(p - q)p and Q'rj := ^(1 + CC)|? - ^(p - q)v, respectively. The 
components of the connection 1-form A e in the NP basis are denoted by 70c, 711', 7oi' and 
710', respectively. The components of the curvature are Xo '■= F a j ) l a m b , 2\\ '■= F a b(l a n b + 
fh a m b ), X2 '■= F a f,rh a n b , Xo '■= F a bl a rh b , 2xi '■= F a b(l a n b + m a m b ) and xi := F a bm a n b . 
Note that for real connections 70c and 711/ are real and 710' is the complex conjugate of 
7oi'. Similarly, also for real connections, Xii X2 and X3 are complex conjugate of Xo, Xi 
and X21 respectively, but for general complex connections they are independent. 

In the gauge 70c = of Newman [3] the expression of the components of the field 
strength by the components of the connection 1-forms and the field equations (12. ip to- 
gether with the Bianchi identities take the form 



dim' 1 ( - 

Xo = ^— + -7oi', (5.1) 
or r 

2X1 = ~&r~ ~ + ^° 3 ' 701 ' + bio',7oi'] , (5.2) 

^T" + -Xi - -o3'xo = - [Xo, 7io'l + 27rjoo', (5.3) 

-fijr + - ^o3'xi = - [xi, 7io'] + 27rj 10 ' , (5.4) 

-g^- = -X2 + 2 -Qf~ + 2^7io' + - 3 7ii' + [lio' ,7ii'J, (5-5) 

"^7 = 2^ £ + 27 Xo+ r° 9Xl+ [ Xo ' 711 '] ~ fci'Toi'] - 2 ^oi', (5.6) 
^ = 2^ + -Xi + -o3X2+ [Xi,7n'] - [X2,70i'] -27Tjn/. (5.7) 

The first four are only the constraint (or 'hypersurface') equations, while the remaining 
three are the evolution equations. In the complex case there are analogous equations for 
the independent curvature components ^0, Xi and X2, too. 



5.1 The solution of the hypersurface equations near a point 

In the small sphere limit calculations of subsection 13.4.11 we need to know only Xi and 
7oi' on the light cone of the point p G M. Thus we need to solve only (I5.ip - fl5.3p with the 
accuracy required by the small sphere calculations. Hence we expand the components 
of the curvature and the connection 1-form as x = X + rx^ + r 2 x^ + 0(r 3 ) and 
j — _)_... _|_ r 3^,(3) _j_ (9( r 4 ^ respectively, as well as j = + rj^ 1 ' + 0(r 2 ). Then the 
solution of (15. ip and its tilded version is 
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7$=0, 7S! = i X r, 1® = ^, lffi = ixf, (5.8) 
7$ = 0, 7S = ixf, 7S = ixi", 7$ = J*?', (M) 
The sum of ( 15. 2 p and its tilded version yield 

1% = XS 0) + XS 0) , 7$ = \(X? + #), 7!? = \(x? + X?)- (5.10) 

Equations ( l5.8p -( l5TTUj) give all the components of the connection 1-form in terms of the 
curvature, except 7^), which remains undetermined. This, however, can be made zero by 
an appropriate gauge transformation. The various components of the curvature are not 
independent either, as they are restricted by (15.31) and (15. 4p . Since we need only \ii we 
evaluate only the former. Using (15. 8p it gives 



XS 0) = "erf, (5.11) 
X? = la#xP + 2 i& (5.12) 
^ 2) = ^ 2) -^Xr,xS 0) ] + |iS, (5.13) 

The difference of these and their tilded version reproduces the expressions for Xo^ ~ Xo\ 
... etc obtained from (15. 2p and its tilded version. The expressions for the expansion 
coefficients of \ 2 = X2 + r %2 1 ' ) + ' ' ' can ^ e derived from (15 .4p in an analogous way. 

5.2 The asymptotics of Yang Mills fields near null infinity 

In [3j Newman solved the source free Yang-Mills equations near the future null infinity, 
and determined the freely specifiable functions of such solutions. In subsection 15.2.11 we 
present the results of the analogous analysis with the higher accuracy of the asymptotic 
expansion required by the analysis of section HI In subsection 15.2.21 the asymptotic field 
equations are solved in certain special stationary configurations. In subsection 15.2.31 we 
discuss some examples with specific gauge groups. 

5.2.1 The asymptotic solution of the Yang— Mills equations 

As Newman showed [3], the hypersurface equations can be integrated in the gauge 700' = 
0, and if we assume that the curvature components xo an d Xo have the asymptotic form 

Xo = ^3X0 + liXl + ^Xl + o{r~ b ) , xo = ^Xo + liXl + ^xl + o(r" 5 ) , 
then the evolution equations imply, in particular, the peeling 

* = ^Xl + + ^Xl + 0(0, X, = \xl + jU + ^Xl + 0(r-% 
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and similar peeling for xi an d X2, too. The functions in these expansions satisfy 



(o37°+ [7°, f}), 
*8 = o»X?+[7°,x;], 

X £ = - 3' ( 3 X ° + [/, X°o] )-[?, oZx°o + [/, X°o}]+3 [xl, x8] , 

= -o3(o3'xS+ [7°,XS] + [xlx°o]) ~ [l°M+[l°,xl] + [x°o,X°o] 
3xS+ [xS,3xS + 2x?] - [xg, 3(^x8 + [7°,x?]) + (o»*8+ [t°,x8]) 



2Xo 



(5.14) 
(5.15) 
(5.16) 

(5.17) 



(5.1* 



and there are the analogous ones for xl, Xi> ■■■ e tc obtained from fl5.14p - fl5.18p by formally 
taking their complex conjugate. (The dot denotes derivative with respect to the retarded 
time coordinate u.) Here 7 = -y°(u, (, C) is an arbitrary function (and, for complex bundle 
and connection, 7 = 7°(it, C> C) is also), but the functions xb — Xo( u iC>C)i i = 0)1)2, 
(and, in the complex case, Xqi too) can be specified freely only at a given retarded time 
instant u = uq. Their time evolution is governed by fl5.16l) - fl5.18p (and the analogous 
tilded equations). Similarly, the sum xl + Xi can be freely specified at u — Uo, but the 
difference xl ~ xl I s already determined by 7 and 7 via 



xl 



xl 



o3 / 7 -o37 +[7°,7 ]. (5.19) 

It is a simple calculation to check that the evolution equation (I5.15P and its tilded version 
preserve f)5.19p . Geometrically 7 and 7 are the leading (r _1 order) term in the asymp- 
totic expansion of 701' and 7 1( y, respectively, and they are analogous to the asymptotic 
shear in general relativity, while xl an d xl represent the outgoing radiation. In terms of 
these the expansion coefficients x\i x\i x\ an d x\ are given by 



x{ = -(oZ'x o+b ,X°o\), (5.20) 

X^-l^'xl+ifiXtt + ixlxt]), (5.21) 

xl = -{oVxl+[l°,xl}), (5.22) 

Xl = ±(o3Vx8+ [7°, 41) + [7° 5 o3'x8+ [t°,x8]] - [*$,*8]), (5.23) 

and there are analogous expressions for x\i Xv X2 an d xhi too. The gauge is not com- 
pletely fixed: gauge transformations with arbitrary A A b = A A b(C> C) can still be carried 
out. 



5.2.2 Special stationary configurations 

Let the Yang-Mills field be stationary with respect to u, i.e. all the dot-derivatives are 
vanishing. Newman and Penrose showed |9] that in the analogous gravitational case the 
asymptotic shear can be made zero by an appropriate supertranslation of the origin cut. 
However, now the integrability condition of the local existence of a gauge transformation 
A(£, C) taking 7 and 7 to zero is the vanishing of the corresponding curvature, which 
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is proportional to the right hand side of (I5.19JI . Hence, for example, in the presence of 
magnetic charges (see equation (I3.2|l ) there is no such a gauge transformation. 

Thus in the rest of this subsection we consider those special stationary solutions of the 
asymptotic field equations in which 7 = and 7 = can be ensured even globally, and 
hence Xi = Xx holds. Thus the total magnetic charges are vanishing. (If the connection 
is real, then Xi is a l so real.) Then equation ( I5.16P immediately gives that Xi is constant 
on S, and hence has the form Xi — \^ e aB f° r some constants C M . These constants are 
real for real bundles and connections, while in the complex case they are complex but 

_ (jv^ gaw j n subsection 13.4.31 that the constants C 1 represent the total charges. 

To solve equation ( 15.1 7p . let us write Xo = h a e^ B for some functions h a with spin 
weight one. Then ( I5.17P takes the simple form 

2 5' 5^ = 3(C a c^)h v =: 3C» u h u . (5.24) 

Thus C^ u is a constant dim g x dim g real or complex matrix. Expanding in terms of 
the s = 1 spin weighted spherical harmonics with the constant coefficients hj m as 

00 j 

k " = E E h "mlYjm (5.25) 
j=l m=-j 

and recalling that the action of the edth operators on spin-s spherical harmonics is 



o3«*Jm = -77= Y {j - s ) (j + s + l)s+lYj m , $sY jm = -7=\ {j + 8 ){j - 8 + l)--l*im, 

equation (I5.24p leads to 

C^h» m = -~(j-l)(j + 2)h>] m . (5.26) 

Thus the expansion coefficients hj m must be eigenvectors of C^ v with the real non-positive 
eigenvalues — |(j — l)(j + 2). By the definition of the charge matrix C^ V C V = holds, i.e. 
A = is always an eigenvalue with eigenvector C M . Thus, in particular, h^ m can always be 
non-zero, and if C 1 is the only eigenvector of C^ v with zero eigenvalue then h^ m = C^h m 
for some constants h m . If the dimension of the kernel space of C fl u is k > 1, then h^ m will 
be the linear combination of the vectors of a basis of this kernel with arbitrary constants 
h m , h m , ... etc. However, since C^ v is a finite matrix, it may have only finitely many 
eigenvalues, and hence only finitely many coefficients hj m can be non-zero. Nevertheless, 
since the total (electric) charges can be specified freely, in general no eigenvalues of 
C^y has the form — — l)(j + 2) for some j > 2. Therefore, apart from the exceptional 
cases in which the charges have this structure, the expansion coefficient functions Xo an d 
Xo are combinations only of the spherical harmonics {¥i m and _iYJ. m , respectively. 

Similarly, if we write Xo = ^ ae aBi then equation (I5.18P reduces to the linear, but 
inhomogeneous equation 

oW + 2F + h -C\k v = -cj, ( ( dh a ) h p + 2h a ( bh p ) + 3h a ( 8'^)) . (5.27) 
Following the strategy of the previous paragraph and expanding k^ as 
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oo j 

k " = E E ( 5 - 28 ) 

j=l m =-j 

for the expansion coefficients kj m of the solution of the homogeneous equation we obtain 
the condition 

C\k] m = 1(0 -I) (J + 2) - 4)*£ m . (5.29) 

Thus the general solution of the homogeneous equation contains non-trivial kj m for j 2 
only if the matrix is chosen to be very special. This solution depends only on finitely 
many parameters. A particular solution of the inhomogeneous equation ( I5.27P can also be 
found by using the spin weighted spherical harmonics and the expansion of their product, 
s Yj m s iYji m ', by the appropriate Clebsch-Gordan coefficients and the spherical harmonics 
s+s rY Jm+m/ , max{|s + s'\, \j - j'\} < J < j + f. 

Thus, to summarize, in the present special stationary case the asymptotic field equa- 
tions can be integrated completely and these solutions (with the second order accuracy, 
i.e. keeping only %q, Xo an d Xi) can be indexed by the freely specifiable total (electric) 
charges and the expansion coefficients h^ m and fc^m ( an d, in the complex case, by 
h^ m and k^, too). Other parameters may occur in the solutions (viz. in the expansion 
coefficients Xo, Xo> e ^ c ) on ^V wnen the charges are chosen in very specific ways. 

5.2.3 Special stationary solutions: Examples 

First consider the generic case, i.e. when no non-vanishing eigenvalue of C M „ has the form 
±|(j — l)(i + 2) for some j G N. Then by (I5.26P we have that = J2m=-i ^lmi^im an d 
= Y2m=-i him-iYim, and then, using the expressions how the edth operators act on 
the spin weighted spherical harmonics, (15. 2 7ft reduces to 

,00 j 

j=l m=-j 
1 

+ Cp E {3hty im ,-2h a lm ht,)iY lm0 Y lm >. (5.30) 

m,m'=— 1 

Its second line can be written as Y,L=-i #imi*im+Em=-2 H-im^m, where H? m , j = 1, 2, 
are explicitly given, homogeneous quadratic expressions of hi m and h^ m) and are built also 
from the structure constants and the Clebsch-Gordan coefficients. With this substitution 
we obtain that 

4 2 2 

(C^y + -5^)^ m = — -Hi m , C^vk^m = — -i^2m- (5.31) 

Moreover, if C^ v has an eigenvalue \{{J — 1)(J + 2) — 4) for some J > 3 then kj m are 
its eigenvectors while kj m = for all j = 3, J — 1, J + 1, ...; otherwise the coefficients 
kj m are vanishing for all j > 3. The first of (I5.3ip can be inverted to find the unique 
solution ki m precisely when — | is not an eigenvalue of C^ V) otherwise the corresponding 
eigenvector can be freely added to the particular solution. Similarly, is determined 
only up to the general solution of the homogeneous equation C^^k^ = 0. 
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If the only eigenvector of C^ v with zero eigenvalue is the charge vector C M , then U{ m = 
C^h m and h^ m = C^h m , which imply that Hj m = 0. (For the sake of simplicity we assume 
that the structure constants are real.) Hence k^ m is either an arbitrary combination of 
the eigenvectors of C^ v with eigenvalue — | or zero, while & 2m is a solution of C^Jz^ = 
0. If, in addition to C M , (7 M is another eigenvector of C^ v with zero eigenvalue, then 
^ = X]m=-i(C^m + C^h m ) iYi m , implying that Hf* m are proportional to cf^ p C a C^ = 
C^pCJP = 0, i.e. Hj m = again. In the other extreme special case when the charge matrix 
is vanishing (e.g. if C M = 0) k^ m = — |iff m holds, k^ m are arbitrary, but now H% m = 
must hold. This yields restrictions on the coefficients h^ m and h^ m . The evaluation of 
the explicit formulae show that then = 0, but H± ±l are still non-zero. 

To illustrate the above possibilities, let us specify the Lie algebra g. First it is chosen 
to be the abstract Lie algebra s/(2,R), so(3) (i.e. su(2)) or s/(2,C); and let us choose 
the basis {ei,e 2 ,e 3 } such that [e 2 ,e 3 ] = T^i, [e 3 , ei] = e 2 and [ei,e 2 ] = e 3 , where the 
upper sign corresponds to s/(2,R) and the lower to so(3) and s/(2,C). In this basis the 
components of the Killing-Cartan metric are g a p = 2diag(— 1, ±1, ±1). The eigenvalues 
of C^ v are A = 0, in which case the corresponding eigenvector is proportional to the total 
charge vector C M , and A = 1 g a ^C a C^ . The non-zero eigenvalues can be real for 

g = s/(2,R) if g a ^C a C^ > and for g = s/(2,C) also in special cases (e.g. for purely 
imaginary C M ), but never for g = so(3). Thus the configurations for g = so(3) are always 
generic. 

From the point of view of general relativity the most interesting choice for g is the 
abstract Lorentz Lie algebra so(l, 3). We choose its basis {ei, • • • , e^} to be the standard 
one, i.e. in which {ei,e 2 ,e3} is the basis above for so(3) while e^, and are the 
'boost generators' such that [e^ej = 0, [ei,es] = e^, ... etc. In this basis the Killing- 
Cartan metric is g a p = 4diag(— 1, —1, — 1, 1, 1, 1). One eigenvalue of the charge matrix is 
zero with eigenvectors and := (— C 4 , —C 5 , — C 6 , C 1 , C 2 , C 3 ) (as a column vector). 
Clearly, g a pC a C^ = —g a pC a C" . The (square of the) non-zero eigenvalues are A 2 = 
\(g a pC a C^ '±i\g a pC a CP\). Thus A is real precisely when g aP C a C p = and g a pC a C p > 0, 
and, in contrast to so(3), for special values of the non-zero eigenvalues could have the 
form ±|(j - l)(j + 2) or - l)(j + 2) - 4) for some j G N. 

Finally we summarize the results in an exceptional case, i.e. when ±|(J — 1)( J + 2) 
are eigenvalues of C^ v for some integer J > 2. Solving (I5.26p . we have that 

l J 

h" = E h lm lYlrn + E h Jm ^ (5-32) 

m=— 1 m=—J 

where hi m and hj m are arbitrary combinations of the eigenvectors of C^ v with zero 
and — g(J — 1)(J + 2) eigenvalue, respectively. If we substitute this and the analogous 
expression for /i M into (I5.27p . use elementary properties of the edth operators and the spin 
weighted spherical harmonics and write the products of such harmonics as combinations 
of spherical harmonics with Clebsch-Gordan coefficients, we obtain 



0= 2^ E ^C^ m +(4-(j-l)(j + 2))kQ 1 Y jm + (5.33) 

j=l m=-j 

1 j 2J j 



m,m'=— 1 



j=l m=-j 
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Here can be given explicitly in terms of h% m , /i" m , hj m , hj m , the structure constants 
and Clebsch-Gordan coefficients. Since we assumed that C^ v has non-zero eigenvalues, it 
cannot be identically zero, and, for the sake of simplicity, let us assume that C^ v does not 
have more than two eigenvectors with zero eigenvalue. Then (%g(3hi m h^ m , — hi m h^ m ,) = 
(see the paragraph following equation (I5.3ip ). Then, evaluating (15.33 p . we find that 
kjm = for j > 2J + 1, kj m are uniquely determined for j = 1, 3, 4, 2 J, and k% m 
are determined up to the addition of the general solution of the homogeneous equation 

C^uk^m — 0. 
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